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The equations of motion for a massive spin-| Rarita-Schwinger field in a 
. finite-range, central, Lorentz scalar potential are developed. It is shown that 
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the resulting density may not be everywhere positive definite. 
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|>. : I. INTRODUCTION 

ON ■ 

^T* \ For many years there has been interest in studying spin-| particles in as much detail as 

has been done for spin-7j particles. In the latter case, the Dirac equation has been shown to 
provide a reasonable starting basis for such studies 0, although of course the formalism of 

• i— i . 

quantum field theory is needed for any detailed consideration. Jasiak and Szymacha @ have 
attempted to use the Rarita-Schwinger equations for spin-| particles in order to investigate 
the A resonance in the bag model. However, they restricted themselves to solving only a 
subset of the equations, and made no attempt to verify whether or not the wave function 
that they calculated was physically meaningful. In this paper, we address this question by 
explicitly considering the most general form for the wave function. 

Attempts to quantize the spin-| Rarita-Schwinger field in the presence of interactions 
of scalar and vector type are beset with problems ||]. For example, it is known that the 
Rarita-Schwinger field has noncausal modes of propagation ||. Johnson and Sudarshan || 
found that for vector interactions, the canonical anti-commutator is not positive definite at 
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all spacetime points. Here we will restrict our considerations to the classical field equations. 

Our work here is motivated by relativistic QHD models of Walecka and collaborators 
|1]U, which consider nucleons interacting with mean- field scalar and vector potentials by 
the Dirac equation. Our initial purpose in undertaking this work was to describe the A 
by an analogous relativistic equation. This is clearly important in order to obtain a full 
treatment of the electromagnetic response of the nucleus, particularly for deep inelastic 
electron scattering in the quasi-elastic regime and beyond f7|. An alternative approach, 
taken by Lin and Serot ||, is to treat the A dynamically as a n — N resonance, thus 
avoiding the need for an explicit equation for the A. Here we treat the A as an additional 
degree of freedom, described by the relativistic spin-| Rarita-Schwinger equation. This is 
analogous to the non-relativistic A-hole models. 

II. FORMALISM FOR NONUNIFORM, CENTRAL SCALAR POTENTIALS 

The free-particle Rarita-Schwinger equations (i ft— m)^ = 0, and the auxiliary condition 
7 M \I/ M = 0, can both be derived from the Lagrangian 

C RS = + i m)tf" - ~* M (<^$ + Yd ■ *) + ^^(P - i m)$, (1) 

where ft = 7^9^, $ = 7 A1 X I /M , and d ■ = d^^. The Euler-Lagrange equations of motion 
corresponding to Crs are given by 

{p + % m)W - + ■ *) + -Y(P - i m)$ = 0. (2) 

3 3 

Contracting this result with 7 M and results in the subsidiary conditions $ = 0, and 
d-ty = 0. Hence Eq. reduces to the usual Rarita-Schwinger equations. Note that Eq. ([[]) 
is not the most general form for £rs- In fact, there is a two-parameter Lagrangian that 
leads to exactly the same free equations of motion regardless of the values of the parameters. 

The most general form of bilinear, non-derivative couplings of a massive Rarita-Schwinger 
field was investigated by Hagen and Singh ||. Here we are interested in the simplest case of 
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a Lorentz scalar potential. In order to introduce this interaction potential into the Rarita- 
Schwinger formalism, we follow the approach used for the Dirac equation (and applied to 
the spin- 1 case by Jasiak and Szymacha [[|) and replace the constant mass m = ttlq in the 
Lagrangian and resulting Euler-Lagrange equations by a mass m(r) = mo — U(r), which is a 
spherically symmetric function of radius r. This approach is consistent with that of Hagen 
and Singh ||. The scalar potential U(r) is of finite range, and vanishes at sufficiently large 
values of r. 

The subsidiary conditions that were derived from Eq. (||) now become 

1 dm . , , 

d-^ + — f * = 0, 3 

m dr 

2d- ^ - zm$ = 0. (4) 

The Euler-Lagrange equations can be written in the forms: 

(^ + im)(f° - - 7 °$) + -(2d° -im7°)$ = 0, (5) 
3 6 

1 9$ 1 1 

r-{p + i m)* + Yd ■ $ + -f(P ~ i rn)& = 0, (6) 

3 ut 3 3 

where Y' = r ■ 7. These two equations, along with the two subsidiary conditions (d) and 
(d), replace the original four Rarita-Schwinger equations. 

Using the usual Dirac representation of the gamma matrices 7 M , it follows that ft = 
7°d° + Y (f: + i±*), where k = -(1 + a ■ L). Then Eq. (§) reduces to 

(Yd + im)(mW)+ (IL + ^LH] (mY^l 

\or r J 

- \ (l°Yd° - ^) (m d>) - ^7° (*° " ^7°$) = 0, (7) 

where \l/ r = f • 

In order to solve these equations, we seek normal-mode solutions for \I/ M (r,i) in their 
most general form consistent with conservation of angular momentum. For the time-like 
component, we have 

( G n (r)iYHr) ^ 



V°(r,t) =e~ iEt - 
r 



K -F {r)Y$(r) j 



(8) 



where Y/(r) = \Ye(r) ® xP, X is the spin-- spinor, and £' = 2j — £ is the other value of £ 
possible for the given value of the angular momentum quantum number j. For the space-like 
component, we have the superposition 



*(r,t) = e" 



-iEt 



G 1 (r) Y{ 3 (r) + G 2 (r) , (*) + G 3 (r) Y{ tl _{r) 
F x (r) i Yj 3 (f) + F 2 (r) i Y) , (f ) + F 3 (r) i (f ) 

£1 2 £2 2 c 2 



\ 
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where Y\ s {r) = \Ye(r) (g> [x <g> ei] s p', with ej the unit vector. The values of £i,£2,£i and £2 
in Eq. (^]) are constrained by angular momentum considerations (£i — j + |), and by the 
Rarita-Schwinger equations themselves, as will be seen shortly. 

Before trying to solve Eqs. (|5[) and ([?]), we first consider several identities for a vector 
field field V{ s {r) = h(r)Y j u (r). These are 



r-Vi(r) 



(10a) 
(10b) 
(10c) 



where 



iff = £+l; 
(£+1), iff = £-l; 



a, 



j+s-if £ 1 £'\f £ I £' 







(11a) 
(lib) 



with [j] = ^2 j + 1 . It then follows from putting Eqs. @ and @ into Eq. ((D that 



$ = e 



-iEi 



{G (r)-V3F 3 (r))iYj 
(F (r) + V3G 3 (r))Yj, 



(12) 



Hence Eq. (|) reduces to the set 



(2E - m)Go(r) - mV3F 3 (r) + |( ^- - £ j ( 2 ^o(r) - >/3G 3 (r)) = 0, 



(13a) 



-(2£ + m)Fo(r) + mV3G 3 (r) + ~ + ~) (2G (r) + ^F 3 (r)) = 0, (13b) 



with k = -(£ + 1) if j = £ + ~, and k = £ if j = £ - \. 

For the other three sets of equations, it is useful to define the following quantities: 



GAt) 



2.7-1 
3 

3(j+l) 



Gi{r) 



Go(r) 



2j+3 
3j 

0+1) 



Go(r) 



FAr) 



2.7-1 
3(i+l) 

3(i+i) 



Fi(r), 



F 2 (r) 



2i+3 
J+l 

H±3 

3j 



Fo(r) 



(14) 



a 



-(i-i: 



a = 



-0' - 1) 



■0" + 1) 

a + 1) 



where the upper (lower) component applies for j = £+\ (j — I — |). We also introduce the 
following four functions: 



Then Eq. 



G(r) = G 1 (r)-G 2 (r)--^G 3 (r), 

F{r) = F 1 (r) - F 2 (r) - ^F 3 {r), 

&(r) = aG^r) - (1 - a)G 2 {r) + A^G 3 (r), 

F'(r) ee a'F 1 (r) - (1 - a')F 2 (r) - 4=Fs(r). 

together with (f|), implies that $ = ^r^^ r , and so we find 

G (r)-V3F 3 (r) = ——G(r), 
m z dr 



F (r) + v / 3G 3 (r) 



2 dm ~ , , 
-F(r) 



m 2 <ir 

The angular momentum constraints for £ = j — - require I' — j ' + ~, l\ — j 



ii=3 



± £ 2 = j + f , £i = j 



and £2 = j + h- Eq. @ reduces to 



^ , n /— , ^ dG(r) G'(r) 

{2E + m G r - mV3F 3 r + 2—^ + 2^^ 

ar r 

(2£7 - m)Fo(r) - m v / 3G 3 (r) + + 2^^ 

ar r 



0, 



0. 



(15a) 

(15b) 
(15c) 
(15d) 

(16a) 
(16b) 



-§.4 = i + 5 
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?i = j — n and £2 = j ' + I. Similarly, the constraints for £ = j + 5 require = j — ^ 
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(17a) 
(17b) 



Finally, Eq. (|7D implies that 

~(E- m)(mG(r)) - ^(2G (r) + V3F 3 (r)) - l -E m(F (r) + V3G 3 (r)) 

+ ^(1 - «)(G? (r) - >/3F s (r)) - f ^ + ^) (mF(r)) = 0, (18a) 
-{E + m)K(r)) - ^(2F (r) - VSG 3 (r)) + ^m(G (r) - V3F 3 (r)) 

+ ^(1 + K )(F (r) + v/3G 3 (r)) + ^ + ^) (mG(r)) = 0. (18b) 

Thus we have derived eight equations (|TB| ^-[T^), in eight unknown functions, that 
together constitute the Rarita-Schwinger equations in a finite, spherically symmetric system. 
Eqs. ( [TBI ) can be used to eliminate G 3 (r) and ^(r), while Eqs. ([IT]) give G"( r ) an d F'(r) 
in terms of G*o(r), F (r), G(r), and F(r). We are then left with four linear, first-order 
differential equations involving only four functions. This system of equations (with energy 
eigenvalue E) can be readily solved using standard numerical techniques. 



III. PROBABILITY DENSITY 

We now turn our attention to the probability density for the system. The conserved 
current J M corresponding to the Lagrangian Crs is given by 

J" = q al W a - - ($"4°^ + + -$ f 7V$. (19) 

Explicitly, we have that d ■ J — 0, if ^ satisfies the Euler-Lagrange equations. Then 
p = J° = po — ^ ' ^ can be considered to define the probability density of the spin-| 
particle, where 

Po = - i($t 7 °^° + ^°t 7 °$) + -$t$ 

3 3 
= ^/l 2 {^o (0 + Ff(r)} + ^l^l 2 (r) + Gl(r)) , (20) 

using the explicit forms for \1/ and The quantity Sl/t • can be evaluated from Eq. @, 
leading to 
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p(r) 



4 7rr 2 



X 



[jf E 



£even 

2/ £ j j 



-rrij 







P c {coae){-l) mi+1/2 
[£' 



\Glir 



+ J2 \G l (r)G k (r)[£ i }[e k } 



i,k=l 



2 ^« 4 



+ F 4 (r)F fc (r) ft] [4] (f | / 



2j £ j j 

\ \ e g 

it 4 £ 
000 



^/s ^ 

000 



e e' £ 
000 



(21) 



= ^ {§G|j(r) - G\{r) - G 2 2 (r) + ^F 2 (r) - F 2 (r) - F 2 2 (r)} + higher-£ terms. (22) 

For j = ~, only the above £ = term is allowed. Similarly, only this term contributes for 
j = I at cos# = 

Explicit numerical solution of the Rarita-Schwinger equations, Eqs. (|T3|) and ([18|), as well 
as the dependent Eqs. ( |i6|) and flTFD , for j = | {I — 0, 1) and j — | (£ = 1, 2) were undertaken 
for several physically plausible potentials. In particular, we examined the solutions for a 
particle of mass 1230 MeV in potentials of nuclear range (a few fm.) and depth (a few 
hundred MeV), corresponding to the A in a nuclear environment. Although reasonable 
values of the energy eigenvalue were found, the wave functions typically caused the term in 
brace brackets in Eq. (p2[) to change sign as r varies. This corresponds to the existence of a 
negative probability density, since only the overall sign of J M lacks any physical significance. 
This indicates either that the interpretation of p(r) as being a probability density is invalid, 
or that the Rarita-Schwinger formalism with a scalar potential is intrinsically unphysical. 
In addition, our attempts to use the two degrees of freedom inherent in the more general 
form for £^5 to see if another choice of the parameters could lead to a p(r) that is explicitly 
positive definite were also unsuccessful. All parameter sets led to the same qualitative 
behavior as was observed for the Lagrangian of Eq. ([!]) . 

The lack of a positive definite probability density suggests that many of the problems 
that plague attempts to quantize spin-| fields are also to be found in the classical field 
equations. 
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